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We show that linearly polarized gluons inside unpolarized hadrons can be directly probed in jet 
or heavy quark pair production in electron-hadron collisions. We discuss the simplest cos2<^> asym- 
metries and estimate their maximal value, concluding that measurements of the unknown linearly 
polarized gluon distribution in the proton should be feasible in future EIC or LHeC experiments. 
Analogous asymmetries in hadron-hadron collisions suffer from factorization breaking contributions 
and would allow to quantify the importance of initial and final state interactions. 
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Although quarks and gluons are confined within 
hadrons, tests of their fundamental properties are pos- 
sible through scattering processes. It has become clear 
that quarks are in general spin-polarized, even within un- 
polarized hadrons, with polarization directions and mag- 
nitudes that depend on their transverse momentum and 
flavor. This nontrivial feature of hadron structure shows 
itself through specific angular asymmetries in scattering 
processes [THHj that have been studied in a number of ex- 
periments |3HZ]- Quark spin-polarization in unpolarized 
hadrons is also supported by first-principle lattice QCD 
calculations [5] . What has received much less attention is 
that gluons can exhibit a similar property, i.e. they can 
be linearly polarized inside an unpolarized hadron. In 
this letter we propose measurements which are directly 
sensitive to this unexplored gluon distribution. Its accu- 
rate measurement would allow one to take advantage of 
polarized scattering at colliders without polarized beams. 

Thus far experimental and theoretical investigations 
of gluons inside hadrons have focussed on their momen- 
tum and helicity distributions. The gluon density g{x) 
describing the distribution of unpolarized gluons with a 
collinear momentum fraction x in an unpolarized hadron, 
integrated over transverse momentum p T , has been ex- 
tracted with considerable precision from measurements of 
high energy electron-proton collisions at HERA (DESY). 
This distribution enters the structure function Ft, in in- 
clusive deep inelastic scattering (DIS) at order a s , and 
it drives the evolution of sea quark distributions at small 
values of x. The unintegrated gluon distribution g(x,p^) 
enters less inclusive reactions where the transverse mo- 
mentum of the gluons is taken into account, such as semi- 
inclusive deep inelastic scattering or dijet production in 
hadronic collisions. In these cases the gluons are not 
necessarily unpolarized, even if the parent hadron itself 
is unpolarized. In fact, because of their spin-orbit cou- 
plings, the gluons can obtain a linear polarization. This 



gives access to a different polarization mode compared to 
the helicity distribution Ag(x), which is the distribution 
of circularly polarized gluons inside polarized nucleons. 

Information on linearly polarized gluons in a hadron 
is formally encoded in the hadron matrix element of a 
correlator of the gluon field strengths F^(0) and F va {\) 
evaluated at fixed light-front time A + = X-n = 0, where n 
is a lightlike vector conjugate to the parent hadron's four- 
momentum P. Specifically, the gluon content of an un- 
polarized hadron at leading twist (omitting gauge links) 
for a gluon momentum p = x P + p T + p~n is described 
by the correlator [5] 

ffi^f, n ) - npn ° [ d ( X - P ) d2X T Ap-X 
9 1 ,Pt) {p-nf J (2tt) 3 

x(P|Tr[F^(0)F^(A)]|F)J LF 

with p\ = -p 2 T , g£ u = g» v - P^n v /P-n - n^P v jP-n. 
This defines the transverse momentum dependent distri- 
bution functions (TMDs) ff{x,p^) representing the un- 
polarized gluon distribution g(x,p^), at fixed light-front 
time, whereas hf 9 (x,p^) is the distribution of linearly 
polarized gluons in an unpolarized hadron. It is named 
h^ 9 , because of its resemblance to the transversely po- 
larized quark distribution inside an unpolarized hadron 
h\ q (also frequently referred to as Boer-Mulders func- 
tion) [J. There are notable differences though: the T- 
odd distribution h ^ q for quarks is a chiral-odd distribu- 
tion (chirality-flip) and it is also odd in p T (it enters as a 
rank 1 tensor). It is zero in the absence of initial or final 
state interactions (ISI/FSI) [ITiHT!?] . The T-even distri- 
bution hf 9 for gluons describes a AL = 2 helicity-flip 
distribution, through a second rank tensor in the relative 
transverse momentum p T (p T -even) . Since an imaginary 
phase is not required for T-even functions, it can in prin- 
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ciple be nonzero in the absence of ISI or FSI. Neverthe- 
less, as any TMD, h^ 9 can receive contributions from ISI 
or FSI, leading to process-dependent gauge links in Eq. 
|lj. Therefore h^ 9 can be non- universal and its extrac- 
tion can be hampered in nonfactorizing cases. 

Thus far no experimental studies of the function h^ 9 
have been performed. It has been pointed out [TS] that it 
contributes to the so-called dijet imbalance in hadronic 
collisions, which is commonly used to extract the average 
partonic intrinsic transverse momentum. Here it enters 
the observable as a convolution of two h^ g functions, 
similarly to the double Boer-Mulders effect which leads 
to a large sin 2 6 cos 2<f> term and the leading-twist viola- 
tion of the Lam- Tung relation in Drell-Yan lepton pair 
production [5J [3] . Although in principle it is possible to 
isolate the contribution from the h^ 9 functions by ap- 
propriate weighting of the planar angular distribution, 
that is likely too difficult to do in practice. Moreover, it 
is unclear whether this weighted observable factorizes to 
begin with, because of factorization breaking effects such 
as discussed in Ref. [13] . 

Given its unique nature, it would be very interesting 
to obtain an extraction of hi 9 in a simple and theoret- 
ically safe manner. This turns out to be possible, since 
unlike hi 9 , it does not need to appear in pairs. In this 
letter we will discuss several new ways to probe the linear 
gluon polarization using observables that involve only a 
single h^ 9 . The processes of interest, semi- inclusive DIS 
to two heavy quarks or to two jets, allow for TMD factor- 
ization and hence a safe extraction. Analogous processes 
in proton-proton collisions run into the problem of fac- 
torization breaking. A difference between the extractions 
will allow to quantify the importance of ISI/FSI. 

We first consider the electroproduction of heavy 
quarks, e(£)+h{P) e{i')+Q(K 1 )+Q(K 2 )+X, where 
the four-momenta of the particles are given within brack- 
ets, and the quark-antiquark pair in the final state is al- 
most back-to-back in the plane perpendicular to the di- 
rection of the exchanged photon and hadron. The calcu- 
lation proceeds along the lines explained in Refs. [HI [15] . 
We obtain for the cross section integrated over the angu- 
lar distribution of the back-scattered electron e(£'): 

dyi dy 2 dy dx B d 2 q T d 2 K± 
(l + yx B ) 
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expression involves the standard DIS variables: 

2 where q is the momentum of the virtual 
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This 

Q 2 = -q 

photon, x B — Q 2 /2P ■ q,y = P ■ q/P ■ £ and s 
(£ + P) 2 = 21 ■ P = 2P ■ q/y = Q 2 /x B y. Further- 
more, we have for the jet momenta K 2 ± = —K 2 ± and 
introduced the rapidities yi for the heavy quark (HQ) 
or jet momenta (along photon-target direction). We de- 
note the heavy (anti)quark mass with Mq. For the par- 
tonic subprocess we have p + q = K\ + K2, implying 



z\ + Z2 = 1, where Zi — P ■ Ki/P ■ q. We introduced 
the sum and difference of the transverse HQ or jet mo- 
menta, K± — (Ki± — K 2 ±)/2 and q T = Ki± + K 2 ± with 
\q T \ <C \K±\. In that situation, we can use the approx- 
imate transverse HQ or jet momenta K\x ~ Kj_ and 
K 2 ± « -K± denoting M 2 ± « M\ ^M^+K 2 ± . The az- 
imuthal angles of q T and K± are denoted by <f> T and <p± , 
respectively. The functions A and B in general depend 
on x B ,y,z{= z 2 ),Q 2 /Ml,M 2 /Ml, and q\. 

The explicit expression for the angular independent 
part A involves only /f . We will focus here on the coef- 
ficient B of the cos2(0 T — <f>±) angular distribution and 
we obtain 
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D = D (z,Q 2 /Ml) = 1 + z(l - z)Q 2 /Ml, a(y) = 2 - 
y(2-y),b(y) = (6-y(6-y))/a(y). 

One observes that the magnitude B of the cos 2</> asym- 
metry, where </> = (f> T — <fi±, is determined by hf 9 and that 
if Q 2 and/or Mq are of the same order as K±, the coeffi- 
cient B is not power suppressed. Since hf 9 is completely 
unknown, we estimate the maximum asymmetry that is 
allowed by the bound: 
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that we derived from the spin density matrix given in 
[5] in the way presented in Ref. [H]. The superscript 
(2) denotes the n — 2 transverse moment. Trans- 
verse moments of TMDs are defined as: f^ n \x) = 
J d 2 p T (p 2 /2M 2 ) f{x,p 2 ) (a suitably chosen regular- 
ization is understood, e.g. as discussed in appendix B 
of Q7]). If we select Q 2 = Mq = K 2 JA, y t = y 2 , the 
asymmetry ratio 



/ d 2 q T q 2 cos 2( 
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is maximally around 13%, which we view as encouraging. 

If one keeps the lepton plane angle <f>^, there are other 
azimuthal dependences such as a cos2(<^ — <fi T ), but its 
bound is at least 6 times smaller than on cos 2{4> T — 4>±). 

The cross section for the process eh — > e' jet jet X 
can be calculated in a similar way. The correspond- 
ing expressions can be obtained from Eqs. ^ and Q 
with Mq = 0. One can then also replace the rapidities 
of the outgoing particles, y±, with the pseudo-rapidities 



3 



rji= — In [tan(!#j)] , 6i being the polar angles of the final 
partons in the virtual photon-hadron cms frame. Note 
that A now also receives a contribution from 7*?-) gq, 
leading to somewhat smaller asymmetries. 

Since the observables involve final-state heavy quarks 
or jets, they require high energy colliders, such as a future 
Electron-Ion Collider (EIC) or the Large Hadron electron 
Collider (LHeC) proposed at CERN. It is essential that 
the individual transverse momenta Kix are reconstructed 
with an accuracy 5K± better than the magnitude of the 
sum of the transverse momenta Kxx + K 2 ± = Qt- Thus 
one has to satisfy 5K± <C \q T \ *C \K±\. 

An analogous asymmetry arises in QED, in the 'tri- 
dents' processes ie(p) — > f/j + /i _ e'(p'orX) or [i~ Z — > 
yrliZ [T8H2T]. This could be described by the distribu- 
tion of linearly polarized photons inside a lepton, pro- 
ton, or atom. QCD adds the twist that for gluons inside 
a hadron, ISI or FSI can considerably modify the result 
depending on the process, for example, in HQ produc- 
tion in hadronic collisions: pp — > Q Q X , which can be 
studied at BNL's Relativistic Heavy Ion Collider (RHIC) 
and CERN's LHC, and pp -> Q Q X at Fcrmilab's Teva- 
tron. Since the description involves two TMDs, breaking 
of TMD factorization becomes a relevant issue, cf. fB] 
and references therein. The cross section for the process 
h 1 (P 1 )+h 2 (P 2 ) -> Q(K 1 )+Q(K 2 )+X can be written in a 
way similar to the hadroproduction of two jets discussed 
in Ref. [13] , in the following form 



da a 2 
d yi dy 2 d 2 K 1± d 2 K 2± = sMj 

x A(q 2 ) + B(q 2 T )q% cos 2(. 
+ C(q\)q% cos4(0 T - 1 



(7) 



Besides q%, the terms A, B and C will depend on other, 
often not explicitly indicated, variables as z, Mq/Mj_ 
and momentum fractions Xx, x 2 obtained from X\j 2 = 
( M xx e ±yi +M 2X e ±y2 ) / y/s . 

In the most naive partonic description the terms A, B, 
and C contain convolutions of TMDs. Schematically, 
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Terms with higher powers in Mq/M]_ are left out. In 
Fig.^Jthe origin of the factor Mq/M]_ in the contribution 
of h 1 9 to B is explained. 

The factorized description in terms of TMDs is prob- 
lematic though. In Ref. [14] it was pointed out that for 
hadron or jet pair production in hadron-hadron scatter- 
ing TMD factorization fails. The ISI/FSI will not allow 
a separation of gauge links into the matrix elements of 



the various TMDs. Only in specific simple cases, such 
as the single Sivers effect, one can find weighted expres- 
sions that do allow a factorized result, but with in gen- 
eral different factors for different diagrams in the partonic 
subprocess [221 [25] . Even if this applies to the present 
case for A and B as well, actually two different func- 
tions h 1 L9 ^ 2 \x) (and fi^\x)) will appear, corresponding 
to gluon operators with the color structures f a t, e fd e and 



cde ■ 



respectively 



This is similar to what 



happens for single transverse spin asymmetries (An) in 
heavy quark production processes P5M5] . Because there 
too two different (/ and d type) gluon correlators arise, 
the single-spin asymmetries in D and D meson produc- 
tion are found to be different. However, in the unpo- 
larized scattering case considered in this letter the situ- 
ation is simpler, since only one operator contributes or 
dominates. In the 7*5 — » QQ subprocess only the ma- 
trix element with the / /-structure appears, while in the 
9 9 Q Q subprocess relevant for hadron-hadron colli- 
sions the d d-structure dominates (the //-contribution is 
suppressed by 1/N 2 ). A side remark on p T broadening 
[3"UH3"2] : because of the two different four-gluon opera- 
tors for fx (x) we expect the broadening Ap 2 in SIDIS, 
(Ap^,)dis = (pt)eA — (Pr)ep! to be different from the one 
in hadron-hadron collisions, (Ap 2 ) hh = (p 2 ) p A — (Pt)pp- 
In case weighting does allow for factorized expres- 
sions, we present here the relevant expressions for B = 
Biq^QQ + (Ml /Ml) B"^ QQ , where 



ggg^QQ 



N 2 { ' [ M\ j 
U q9 (xx,x 2 ,q 2 T ) + U 99 (xx,x 2 ,q 2 T ) 

A B 1 U"(xx,x 2 ,q 2 T ), 



N 2 — I 



(8) 





FIG. 1: Examples of subprocesses contributing to the cos 20 
asymmetries in ep — > e'QQX and pp — > QQX, respec- 
tively. As the helicities of the photons and gluons indicate, 
the latter process requires helicity flip in quark propagators 
resulting in an Ml /Ml factor. 
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with TV being the number of quark colors and 
B x = z(l-z) (z 2 + (l-z) 2 - 




Ml 



(9) 



The weighted integrals which appear in the q\jM 2 - 
weighted cross section (cf. [13]), for Mi = M 2 = M are: 



tt / dq\ [j^] q 2 T H™(x 1 ,x 2 ,q 2 T ) = 



lT \M 
8 J2 hi^\ Xl )ht^\x 2 ) 



(10) 



flavors 

already discussed in [13] . and 



* I dql {w) £ n "^ 

4(hi 9(2 \ Xl ) ff(x 2 ) + fS( Xl )ht g(2) (x 2 )) . (11) 



Whether gg — > QQ is more important than qq — 
depends strongly on the values of Xi and Mq/M\, and on 
whether one deals with pp or pp. In pp collisions and for 
K\ not too large compared to Mq , the contribution from 
9 is expected to be the dominant one. The importance 
of the contribution from q can be assessed through a 
comparison to photon-jet production [IS] . 

In summary, measurements of the azimuthal asymme- 
try of jet or heavy quark pair production in ep and in pp 
or pp collisions (and possibly also in diphoton or even 
Higgs production [331 IS]) can directly probe the dis- 
tribution of linearly polarized gluons inside unpolarized 
hadrons. From a theoretical viewpoint this asymmetry in 
the ep process is among the simplest TMD observables. 
Breaking of TMD factorization is expected in the pp case 
and may lead to uncontrolled corrections. A comparison 
between extractions from these two types of processes 
would therefore be very interesting. The relative sim- 
plicity of the proposed measurements (polarized beams 
are not required) suggests a promising prospect for the 
extraction of this gluon distribution in the future and for 
the study of its potential process dependence. 

We thank John Collins, Markus Diehl, Francesco Mur- 
gia, Ted Rogers, and Jianwei Qiu for useful discussions. 
CP. is supported by Regione Autonoma della Sardegna 
(RAS) through a research grant under the PO Sardegna 
FSE 2007-2013, L.R. 7/2007, "Promozione della ricerca 
scientifica e dell'innovazione tecnologica in Sardegna". 
This research is part of the FP7 EU-programme Hadron 
Physics (No. 227431). 



t Electronic address: s.ibth@slac.stanford.edu 

* Electronic address: mulders@few.vu.nl 

5 Electronic address: cristian.pisano@ca.infn.it 

[1] D. Boer and RJ. Mulders, Phys. Rev. D 57, 5780 (1998). 

[2] D. Boer, Phys. Rev. D 60, 014012 (1999). 

[3] D. Boer, S.J. Brodsky and D.S. Hwang, Phys. Rev. D 
67, 054003 (2003). 

[4] S. Falciano et al, Z. Phys. C 31, 513 (1986). 

[5] M. Guanziroli et al, Z. Phys. C 37, 545 (1988). 

[6] J.S. Conway et al, Phys. Rev. D 39, 92 (1989). 

[7] L.Y. Zhu et al, Phys. Rev. Lett. 99, 082301 (2007); Phys. 
Rev. Lett. 102, 182001 (2009). 

[8] D. Brommel et al, Phys. Rev. Lett. 101, 122001 (2008). 

[9] P.J. Mulders and J. Rodrigues, Phys. Rev. D 63, 094021 
(2001). 

[10] J.C. Collins, Phys. Lett. B 536, 43 (2002). 

[11] S.J. Brodsky, D.S. Hwang and I. Schmidt, Phys. Lett. B 

530, 99 (2002); Nucl. Phys. B 642, 344 (2002). 
[12] A.V. Belitsky, X. Ji and F. Yuan, Nucl. Phys. B 656, 165 

(2003). 

[13] D. Boer, P.J. Mulders and C. Pisano, Phys. Rev. D. 80, 
094017 (2009). 

[14] T.C. Rogers and P.J. Mulders, Phys. Rev. D 81, 094006 

(2010) . 

[15] D. Boer, P.J. Mulders and C. Pisano, Phys. Lett. B 660, 
360 (2008). 

[16] A. Bacchetta et al, Phys. Rev. Lett. 85, 712 (2000). 
[17] A. Bacchetta et al, JHEP 0808, 023 (2008). 
[18] J.D. Bjorken and M.C. Chen, Phys. Rev. 154, 1335 
(1966). 

[19] S.J. Brodsky and S.C.C. Ting, Phys. Rev. 145, 1018 
(1966). 

[20] M.J. Tannenbaum, Phys. Rev. 167, 1308 (1968). 

[21] M. Chick, C. Pisano, E. Reya and I. Schienbein, Eur. 

Phys. J. C 27, 427 (2003). 
[22] C.J. Bomhof, P.J. Mulders and F. Pijlman, Eur. Phys. J. 

C 47, 147 (2006). 
[23] C.J. Bomhof and P.J. Mulders, Nucl. Phys. B 795, 409 

(2008). 

[24] F. Dominguez, B-W. Xiao, and F. Yuan, Phys. Rev. Lett. 

106, 022301 (2011). 
[25] F. Yuan, Phys. Rev. D 78, 014024 (2008). 
[26] Z.B. Kang and J. W. Qiu, Phys. Rev. D 78, 034005 

(2008). 

[27] F. Yuan and J. Zhou, Phys. Lett. B 668, 216 (2008). 
[28] Z.B. Kang et al, Phys. Rev. D 78, 114013 (2008). 
[29] M. Anselmino et al, Phys. Rev. D 70, 074025 (2004). 
[30] M. Luo, J. Qiu and G. Sterman, Phys. Rev. D 49, 4493 
(1994). 

[31] X. Guo, Phys. Rev. D 58, 114033 (1998). 
[32] Z.t. Liang, X.N. Wang and J. Zhou, Phys. Rev. D 77, 
125010 (2008). 

[33] P.M. Nadolsky et al, Phys. Rev. D 76, 013008 (2007). 
[34] S. Catani and M. Grazzini, Nucl. Phys. B 845, 297 

(2011) . 



* Electronic address: D.Boer@rug.nl 



